Trajectories of clinical characteristics, complications and treatment choices in data-driven subgroups of type 2 diabetes.

AIMS/HYPOTHESIS: This study aimed to explore the added value of subgroups that categorise individuals with type 2 diabetes by k-means clustering for two primary care registries (the Netherlands and Scotland), inspired by Ahlqvist's novel diabetes subgroups and previously analysed by Slieker et al. METHODS: We used two Dutch and Scottish diabetes cohorts (N=3054 and 6145; median follow-up=11.2 and 12.3 years, respectively) and defined five subgroups by k-means clustering with age at baseline, BMI, HbA1c, HDL-cholesterol and C-peptide. We investigated differences between subgroups by trajectories of risk factor values (random intercept models), time to diabetes-related complications (logrank tests and Cox models) and medication patterns (multinomial logistic models). We also compared directly using the clustering indicators as predictors of progression vs the k-means discrete subgroups. Cluster consistency over follow-up was assessed. RESULTS: Subgroups' risk factors were significantly different, and these differences remained generally consistent over follow-up. Among all subgroups, individuals with severe insulin resistance faced a significantly higher risk of myocardial infarction both before (HR 1.65; 95% CI 1.40, 1.94) and after adjusting for age effect (HR 1.72; 95% CI 1.46, 2.02) compared with mild diabetes with high HDL-cholesterol. Individuals with severe insulin-deficient diabetes were most intensively treated, with more than 25% prescribed insulin at 10 years of diagnosis. For severe insulin-deficient diabetes relative to mild diabetes, the relative risks for using insulin relative to no common treatment would be expected to increase by a factor of 3.07 (95% CI 2.73, 3.44), holding other factors constant. Clustering indicators were better predictors of progression variation relative to subgroups, but prediction accuracy may improve after combining both. Clusters were consistent over 8 years with an accuracy ranging from 59% to 72%. CONCLUSIONS/INTERPRETATION: Data-driven subgroup allocations were generally consistent over follow-up and captured significant differences in risk factor trajectories, medication patterns and complication risks. Subgroups serve better as a complement rather than as a basis for compressing clustering indicators.

Genomics Research of Lifetime Depression in the Netherlands: The BIObanks Netherlands Internet Collaboration (BIONIC) Project.

In this cohort profile article we describe the lifetime major depressive disorder (MDD) database that has been established as part of the BIObanks Netherlands Internet Collaboration (BIONIC). Across the Netherlands we collected data on Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) lifetime MDD diagnosis in 132,850 Dutch individuals. Currently, N = 66,684 of these also have genomewide single nucleotide polymorphism (SNP) data. We initiated this project because the complex genetic basis of MDD requires large population-wide studies with uniform in-depth phenotyping. For standardized phenotyping we developed the LIDAS (LIfetime Depression Assessment Survey), which then was used to measure MDD in 11 Dutch cohorts. Data from these cohorts were combined with diagnostic interview depression data from 5 clinical cohorts to create a dataset of N = 29,650 lifetime MDD cases (22%) meeting DSM-5 criteria and 94,300 screened controls. In addition, genomewide genotype data from the cohorts were assembled into a genomewide association study (GWAS) dataset of N = 66,684 Dutch individuals (25.3% cases). Phenotype data include DSM-5-based MDD diagnoses, sociodemographic variables, information on lifestyle and BMI, characteristics of depressive symptoms and episodes, and psychiatric diagnosis and treatment history. We describe the establishment and harmonization of the BIONIC phenotype and GWAS datasets and provide an overview of the available information and sample characteristics. Our next step is the GWAS of lifetime MDD in the Netherlands, with future plans including fine-grained genetic analyses of depression characteristics, international collaborations and multi-omics studies.

An omics-based machine learning approach to predict diabetes progression: a RHAPSODY study.

AIMS/HYPOTHESIS: People with type 2 diabetes are heterogeneous in their disease trajectory, with some progressing more quickly to insulin initiation than others. Although classical biomarkers such as age, HbA1c and diabetes duration are associated with glycaemic progression, it is unclear how well such variables predict insulin initiation or requirement and whether newly identified markers have added predictive value. METHODS: In two prospective cohort studies as part of IMI-RHAPSODY, we investigated whether clinical variables and three types of molecular markers (metabolites, lipids, proteins) can predict time to insulin requirement using different machine learning approaches (lasso, ridge, GRridge, random forest). Clinical variables included age, sex, HbA1c, HDL-cholesterol and C-peptide. Models were run with unpenalised clinical variables (i.e. always included in the model without weights) or penalised clinical variables, or without clinical variables. Model development was performed in one cohort and the model was applied in a second cohort. Model performance was evaluated using Harrel's C statistic. RESULTS: Of the 585 individuals from the Hoorn Diabetes Care System (DCS) cohort, 69 required insulin during follow-up (1.0-11.4 years); of the 571 individuals in the Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS) cohort, 175 required insulin during follow-up (0.3-11.8 years). Overall, the clinical variables and proteins were selected in the different models most often, followed by the metabolites. The most frequently selected clinical variables were HbA1c (18 of the 36 models, 50%), age (15 models, 41.2%) and C-peptide (15 models, 41.2%). Base models (age, sex, BMI, HbA1c) including only clinical variables performed moderately in both the DCS discovery cohort (C statistic 0.71 [95% CI 0.64, 0.79]) and the GoDARTS replication cohort (C 0.71 [95% CI 0.69, 0.75]). A more extensive model including HDL-cholesterol and C-peptide performed better in both cohorts (DCS, C 0.74 [95% CI 0.67, 0.81]; GoDARTS, C 0.73 [95% CI 0.69, 0.77]). Two proteins, lactadherin and proto-oncogene tyrosine-protein kinase receptor, were most consistently selected and slightly improved model performance. CONCLUSIONS/INTERPRETATION: Using machine learning approaches, we show that insulin requirement risk can be modestly well predicted by predominantly clinical variables. Inclusion of molecular markers improves the prognostic performance beyond that of clinical variables by up to 5%. Such prognostic models could be useful for identifying people with diabetes at high risk of progressing quickly to treatment intensification. DATA AVAILABILITY: Summary statistics of lipidomic, proteomic and metabolomic data are available from a Shiny dashboard at https://rhapdata-app.vital-it.ch .

Circulating metabolomic markers linking diabetic kidney disease and incident cardiovascular disease in type 2 diabetes: analyses from the Hong Kong Diabetes Biobank.

AIMS/HYPOTHESIS: The aim of this study was to describe the metabolome in diabetic kidney disease (DKD) and its association with incident CVD in type 2 diabetes, and identify prognostic biomarkers. METHODS: From a prospective cohort of individuals with type 2 diabetes, baseline sera (N=1991) were quantified for 170 metabolites using NMR spectroscopy with median 5.2 years of follow-up. Associations of chronic kidney disease (CKD, eGFR<60 ml/min per 1.73 m2) or severely increased albuminuria with each metabolite were examined using linear regression, adjusted for confounders and multiplicity. Associations between DKD (CKD or severely increased albuminuria)-related metabolites and incident CVD were examined using Cox regressions. Metabolomic biomarkers were identified and assessed for CVD prediction and replicated in two independent cohorts. RESULTS: At false discovery rate (FDR)<0.05, 156 metabolites were associated with DKD (151 for CKD and 128 for severely increased albuminuria), including apolipoprotein B-containing lipoproteins, HDL, fatty acids, phenylalanine, tyrosine, albumin and glycoprotein acetyls. Over 5.2 years of follow-up, 75 metabolites were associated with incident CVD at FDR<0.05. A model comprising age, sex and three metabolites (albumin, triglycerides in large HDL and phospholipids in small LDL) performed comparably to conventional risk factors (C statistic 0.765 vs 0.762, p=0.893) and adding the three metabolites further improved CVD prediction (C statistic from 0.762 to 0.797, p=0.014) and improved discrimination and reclassification. The 3-metabolite score was validated in independent Chinese and Dutch cohorts. CONCLUSIONS/INTERPRETATION: Altered metabolomic signatures in DKD are associated with incident CVD and improve CVD risk stratification.

Genome-wide association study of serum magnesium in type 2 diabetes.

People with type 2 diabetes have a tenfold higher prevalence of hypomagnesemia, which is suggested to be caused by low dietary magnesium intake, medication use, and genetics. This study aims to identify the genetic loci that influence serum magnesium concentration in 3466 people with type 2 diabetes. The GWAS models were adjusted for age, sex, eGFR, and HbA1c. Associated traits were identified using publicly available data from GTEx consortium, a human kidney eQTL atlas, and the Open GWAS database. The GWAS identified a genome-wide significant locus in TAF3 (p = 2.9 × 10-9) in people with type 2 diabetes. In skeletal muscle, loci located in TAF3 demonstrate an eQTL link to ATP5F1C, a gene that is involved in the formation of Mg2+-ATP. Serum Mg2+ levels were associated with MUC1/TRIM46 (p = 2.9 × 10-7), SHROOM3 (p = 4.0 × 10-7), and SLC22A7 (p = 1.0 × 10-6) at nominal significance, which is in combination with the eQTL data suggesting that they are possible candidates for renal failure. Several genetic loci were in agreement with previous genomic studies which identified MUC1/TRIM46 (Pmeta = 6.9 × 10-29, PQ = 0.81) and SHROOM3 (Pmeta = 2.9 × 10-27, PQ = 0.04) to be associated with serum Mg2+ in the general population. In conclusion, serum magnesium concentrations are associated with genetic variability around the regions of TAF3, MUC1/TRIM46, SHROOM3, and SLC22A7 in type 2 diabetes.

Apolipoprotein-CIII O-Glycosylation, a Link between GALNT2 and Plasma Lipids.

Apolipoprotein-CIII (apo-CIII) is involved in triglyceride-rich lipoprotein metabolism and linked to beta-cell damage, insulin resistance, and cardiovascular disease. Apo-CIII exists in four main proteoforms: non-glycosylated (apo-CIII0a), and glycosylated apo-CIII with zero, one, or two sialic acids (apo-CIII0c, apo-CIII1 and apo-CIII2). Our objective is to determine how apo-CIII glycosylation affects lipid traits and type 2 diabetes prevalence, and to investigate the genetic basis of these relations with a genome-wide association study (GWAS) on apo-CIII glycosylation. We conducted GWAS on the four apo-CIII proteoforms in the DiaGene study in people with and without type 2 diabetes (n = 2318). We investigated the relations of the identified genetic loci and apo-CIII glycosylation with lipids and type 2 diabetes. The associations of the genetic variants with lipids were replicated in the Diabetes Care System (n = 5409). Rs4846913-A, in the GALNT2-gene, was associated with decreased apo-CIII0a. This variant was associated with increased high-density lipoprotein cholesterol and decreased triglycerides, while high apo-CIII0a was associated with raised high-density lipoprotein-cholesterol and triglycerides. Rs67086575-G, located in the IFT172-gene, was associated with decreased apo-CIII2 and with hypertriglyceridemia. In line, apo-CIII2 was associated with low triglycerides. On a genome-wide scale, we confirmed that the GALNT2-gene plays a major role i O-glycosylation of apolipoprotein-CIII, with subsequent associations with lipid parameters. We newly identified the IFT172/NRBP1 region, in the literature previously associated with hypertriglyceridemia, as involved in apolipoprotein-CIII sialylation and hypertriglyceridemia. These results link genomics, glycosylation, and lipid metabolism, and represent a key step towards unravelling the importance of O-glycosylation in health and disease.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits.

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

IgG N-glycans are associated with prevalent and incident complications of type 2 diabetes.

AIMS/HYPOTHESIS: Inflammation is important in the development of type 2 diabetes complications. The N-glycosylation of IgG influences its role in inflammation. To date, the association of plasma IgG N-glycosylation with type 2 diabetes complications has not been extensively investigated. We hypothesised that N-glycosylation of IgG may be related to the development of complications of type 2 diabetes. METHODS: In three independent type 2 diabetes cohorts, plasma IgG N-glycosylation was measured using ultra performance liquid chromatography (DiaGene n = 1815, GenodiabMar n = 640) and mass spectrometry (Hoorn Diabetes Care Study n = 1266). We investigated the associations of IgG N-glycosylation (fucosylation, galactosylation, sialylation and bisection) with incident and prevalent nephropathy, retinopathy and macrovascular disease using Cox- and logistic regression, followed by meta-analyses. The models were adjusted for age and sex and additionally for clinical risk factors. RESULTS: IgG galactosylation was negatively associated with prevalent and incident nephropathy and macrovascular disease after adjustment for clinical risk factors. Sialylation was negatively associated with incident diabetic nephropathy after adjustment for clinical risk factors. For incident retinopathy, similar associations were found for galactosylation, adjusted for age and sex. CONCLUSIONS: We showed that IgG N-glycosylation, particularly galactosylation and to a lesser extent sialylation, is associated with a higher prevalence and future development of macro- and microvascular complications of diabetes. These findings indicate the predictive potential of IgG N-glycosylation in diabetes complications and should be analysed further in additional large cohorts to obtain the power to solidify these conclusions.

Potential Value of Identifying Type 2 Diabetes Subgroups for Guiding Intensive Treatment: A Comparison of Novel Data-Driven Clustering With Risk-Driven Subgroups.

OBJECTIVE: To estimate the impact on lifetime health and economic outcomes of different methods of stratifying individuals with type 2 diabetes, followed by guideline-based treatment intensification targeting BMI and LDL in addition to HbA1c. RESEARCH DESIGN AND METHODS: We divided 2,935 newly diagnosed individuals from the Hoorn Diabetes Care System (DCS) cohort into five Risk Assessment and Progression of Diabetes (RHAPSODY) data-driven clustering subgroups (based on age, BMI, HbA1c, C-peptide, and HDL) and four risk-driven subgroups by using fixed cutoffs for HbA1c and risk of cardiovascular disease based on guidelines. The UK Prospective Diabetes Study Outcomes Model 2 estimated discounted expected lifetime complication costs and quality-adjusted life-years (QALYs) for each subgroup and across all individuals. Gains from treatment intensification were compared with care as usual as observed in DCS. A sensitivity analysis was conducted based on Ahlqvist subgroups. RESULTS: Under care as usual, prognosis in the RHAPSODY data-driven subgroups ranged from 7.9 to 12.6 QALYs. Prognosis in the risk-driven subgroups ranged from 6.8 to 12.0 QALYs. Compared with homogenous type 2 diabetes, treatment for individuals in the high-risk subgroups could cost 22.0% and 25.3% more and still be cost effective for data-driven and risk-driven subgroups, respectively. Targeting BMI and LDL in addition to HbA1c might deliver up to 10-fold increases in QALYs gained. CONCLUSIONS: Risk-driven subgroups better discriminated prognosis. Both stratification methods supported stratified treatment intensification, with the risk-driven subgroups being somewhat better in identifying individuals with the most potential to benefit from intensive treatment. Irrespective of stratification approach, better cholesterol and weight control showed substantial potential for health gains.

Altered blood gene expression in the obesity-related type 2 diabetes cluster may be causally involved in lipid metabolism: a Mendelian randomisation study.

AIMS/HYPOTHESIS: The aim of this study was to identify differentially expressed long non-coding RNAs (lncRNAs) and mRNAs in whole blood of people with type 2 diabetes across five different clusters: severe insulin-deficient diabetes (SIDD), severe insulin-resistant diabetes (SIRD), mild obesity-related diabetes (MOD), mild diabetes (MD) and mild diabetes with high HDL-cholesterol (MDH). This was to increase our understanding of different molecular mechanisms underlying the five putative clusters of type 2 diabetes. METHODS: Participants in the Hoorn Diabetes Care System (DCS) cohort were clustered based on age, BMI, HbA1c, C-peptide and HDL-cholesterol. Whole blood RNA-seq was used to identify differentially expressed lncRNAs and mRNAs in a cluster compared with all others. Differentially expressed genes were validated in the Innovative Medicines Initiative DIabetes REsearCh on patient straTification (IMI DIRECT) study. Expression quantitative trait loci (eQTLs) for differentially expressed RNAs were obtained from a publicly available dataset. To estimate the causal effects of RNAs on traits, a two-sample Mendelian randomisation analysis was performed using public genome-wide association study (GWAS) data. RESULTS: Eleven lncRNAs and 175 mRNAs were differentially expressed in the MOD cluster, the lncRNA AL354696.2 was upregulated in the SIDD cluster and GPR15 mRNA was downregulated in the MDH cluster. mRNAs and lncRNAs that were differentially expressed in the MOD cluster were correlated among each other. Six lncRNAs and 120 mRNAs validated in the IMI DIRECT study. Using two-sample Mendelian randomisation, we found 52 mRNAs to have a causal effect on anthropometric traits (n=23) and lipid metabolism traits (n=10). GPR146 showed a causal effect on plasma HDL-cholesterol levels (p = 2×10-15), without evidence for reverse causality. CONCLUSIONS/INTERPRETATION: Multiple lncRNAs and mRNAs were found to be differentially expressed among clusters and particularly in the MOD cluster. mRNAs in the MOD cluster showed a possible causal effect on anthropometric traits, lipid metabolism traits and blood cell fractions. Together, our results show that individuals in the MOD cluster show aberrant RNA expression of genes that have a suggested causal role on multiple diabetes-relevant traits.

Discovery of drug-omics associations in type 2 diabetes with generative deep-learning models.

The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities.

Assessment of the bi-directional relationship between blood mitochondrial DNA copy number and type 2 diabetes mellitus: a multivariable-adjusted regression and Mendelian randomisation study.

AIMS/HYPOTHESIS: Mitochondrial dysfunction, which can be approximated by blood mitochondrial DNA copy number (mtDNA-CN), has been implicated in the pathogenesis of type 2 diabetes mellitus. Thus far, however, insights from prospective cohort studies and Mendelian randomisation (MR) analyses on this relationship are limited. We assessed the association between blood mtDNA-CN and incident type 2 diabetes using multivariable-adjusted regression analyses, and the associations between blood mtDNA-CN and type 2 diabetes and BMI using bi-directional MR. METHODS: Multivariable-adjusted Cox proportional hazard models were used to estimate the association between blood mtDNA-CN and incident type 2 diabetes in 285,967 unrelated European individuals from UK Biobank free of type 2 diabetes at baseline. Additionally, a cross-sectional analysis was performed to investigate the association between blood mtDNA-CN and BMI. We also assessed the potentially causal relationship between blood mtDNA-CN and type 2 diabetes (N=898,130 from DIAGRAM, N=215,654 from FinnGen) and BMI (N=681,275 from GIANT) using bi-directional two-sample MR. RESULTS: During a median follow-up of 11.87 years, 15,111 participants developed type 2 diabetes. Participants with a higher level of blood mtDNA-CN are at lower risk of developing type 2 diabetes (HR 0.90 [95% CI 0.89, 0.92]). After additional adjustment for BMI and other confounders, these results attenuated moderately and remained present. The multivariable-adjusted cross-sectional analyses showed that higher blood mtDNA-CN was associated with lower BMI (-0.12 [95% CI -0.14, -0.10]) kg/m2. In the bi-directional MR analyses, we found no evidence for causal associations between blood mtDNA-CN and type 2 diabetes, and blood mtDNA-CN and BMI in either direction. CONCLUSIONS/INTERPRETATION: The results from the present study indicate that the observed association between low blood mtDNA-CN and higher risk of type 2 diabetes is likely not causal.

Diabetes risk loci-associated pathways are shared across metabolic tissues.

AIMS/HYPOTHESIS: Numerous genome-wide association studies have been performed to understand the influence of genetic variation on type 2 diabetes etiology. Many identified risk variants are located in non-coding and intergenic regions, which complicates understanding of how genes and their downstream pathways are influenced. An integrative data approach will help to understand the mechanism and consequences of identified risk variants. METHODS: In the current study we use our previously developed method CONQUER to overlap 403 type 2 diabetes risk variants with regulatory, expression and protein data to identify tissue-shared disease-relevant mechanisms. RESULTS: One SNP rs474513 was found to be an expression-, protein- and metabolite QTL. Rs474513 influenced LPA mRNA and protein levels in the pancreas and plasma, respectively. On the pathway level, in investigated tissues most SNPs linked to metabolism. However, in eleven of the twelve tissues investigated nine SNPs were linked to differential expression of the ribosome pathway. Furthermore, seven SNPs were linked to altered expression of genes linked to the immune system. Among them, rs601945 was found to influence multiple HLA genes, including HLA-DQA2, in all twelve tissues investigated. CONCLUSION: Our results show that in addition to the classical metabolism pathways, other pathways may be important to type 2 diabetes that show a potential overlap with type 1 diabetes.

Four groups of type 2 diabetes contribute to the etiological and clinical heterogeneity in newly diagnosed individuals: An IMI DIRECT study.

The presentation and underlying pathophysiology of type 2 diabetes (T2D) is complex and heterogeneous. Recent studies attempted to stratify T2D into distinct subgroups using data-driven approaches, but their clinical utility may be limited if categorical representations of complex phenotypes are suboptimal. We apply a soft-clustering (archetype) method to characterize newly diagnosed T2D based on 32 clinical variables. We assign quantitative clustering scores for individuals and investigate the associations with glycemic deterioration, genetic risk scores, circulating omics biomarkers, and phenotypic stability over 36 months. Four archetype profiles represent dysfunction patterns across combinations of T2D etiological processes and correlate with multiple circulating biomarkers. One archetype associated with obesity, insulin resistance, dyslipidemia, and impaired ß cell glucose sensitivity corresponds with the fastest disease progression and highest demand for anti-diabetic treatment. We demonstrate that clinical heterogeneity in T2D can be mapped to heterogeneity in individual etiological processes, providing a potential route to personalized treatments.

Genome-Wide Meta-analysis Identifies Genetic Variants Associated With Glycemic Response to Sulfonylureas.

OBJECTIVE: Sulfonylureas, the first available drugs for the management of type 2 diabetes, remain widely prescribed today. However, there exists significant variability in glycemic response to treatment. We aimed to establish heritability of sulfonylurea response and identify genetic variants and interacting treatments associated with HbA1c reduction. RESEARCH DESIGN AND METHODS: As an initiative of the Metformin Genetics Plus Consortium (MetGen Plus) and the DIabetes REsearCh on patient straTification (DIRECT) consortium, 5,485 White Europeans with type 2 diabetes treated with sulfonylureas were recruited from six referral centers in Europe and North America. We first estimated heritability using the generalized restricted maximum likelihood approach and then undertook genome-wide association studies of glycemic response to sulfonylureas measured as HbA1c reduction after 12 months of therapy followed by meta-analysis. These results were supported by acute glipizide challenge in humans who were naïve to type 2 diabetes medications, cis expression quantitative trait loci (eQTL), and functional validation in cellular models. Finally, we examined for possible drug-drug-gene interactions. RESULTS: After establishing that sulfonylurea response is heritable (mean ± SEM 37 ± 11%), we identified two independent loci near the GXYLT1 and SLCO1B1 genes associated with HbA1c reduction at a genome-wide scale (P < 5 × 10-8). The C allele at rs1234032, near GXYLT1, was associated with 0.14% (1.5 mmol/mol), P = 2.39 × 10-8), lower reduction in HbA1c. Similarly, the C allele was associated with higher glucose trough levels (ß = 1.61, P = 0.005) in healthy volunteers in the SUGAR-MGH given glipizide (N = 857). In 3,029 human whole blood samples, the C allele is a cis eQTL for increased expression of GXYLT1 (ß = 0.21, P = 2.04 × 10-58). The C allele of rs10770791, in an intronic region of SLCO1B1, was associated with 0.11% (1.2 mmol/mol) greater reduction in HbA1c (P = 4.80 × 10-8). In 1,183 human liver samples, the C allele at rs10770791 is a cis eQTL for reduced SLCO1B1 expression (P = 1.61 × 10-7), which, together with functional studies in cells expressing SLCO1B1, supports a key role for hepatic SLCO1B1 (encoding OATP1B1) in regulation of sulfonylurea transport. Further, a significant interaction between statin use and SLCO1B1 genotype was observed (P = 0.001). In statin nonusers, C allele homozygotes at rs10770791 had a large absolute reduction in HbA1c (0.48 ± 0.12% [5.2 ± 1.26 mmol/mol]), equivalent to that associated with initiation of a dipeptidyl peptidase 4 inhibitor. CONCLUSIONS: We have identified clinically important genetic effects at genome-wide levels of significance, and important drug-drug-gene interactions, which include commonly prescribed statins. With increasing availability of genetic data embedded in clinical records these findings will be important in prescribing glucose-lowering drugs.

Plasma protein N-glycosylation is associated with cardiovascular disease, nephropathy, and retinopathy in type 2 diabetes.

INTRODUCTION: Although associations of total plasma N-glycome (TPNG) with type 2 diabetes have been reported, little is known on the role of TPNG in type 2 diabetes complications, a major cause of type 2 diabetes-related morbidity and mortality. Here, we assessed TPNG in relation to type 2 diabetes complications in subsamples of two Dutch cohorts using mass spectrometry (n=1815 in DiaGene and n=1518 in Hoorn Diabetes Care System). RESEARCH DESIGN AND METHODS: Blood plasma samples and technical replicates were pipetted into 96-well plates in a randomized manner. Peptide:N-glycosidase F (PNGase F) was used to release N-glycans, whereafter sialic acids were derivatized for stabilization and linkage differentiation. After total area normalization, 68 individual glycan compositions were quantified in total and were used to calculate 45 derived traits which reflect structural features of glycosylation. Associations of glycan features with prevalent and incident microvascular or macrovascular complications were tested in logistic and Cox regression in both independent cohorts and the results were meta-analyzed. RESULTS: Our results demonstrated similarities between incident and prevalent complications. The strongest association for prevalent cardiovascular disease was a high level of bisection on a group of diantennary glycans (A2FS0B; OR=1.38, p=1.34×10-11), while for prevalent nephropathy the increase in 2,6-sialylation on triantennary glycans was most pronounced (A3E; OR=1.28, p=9.70×10-6). Several other TPNG features, including fucosylation, galactosylation, and sialylation, firmly demonstrated associations with prevalent and incident complications of type 2 diabetes. CONCLUSIONS: These findings may provide a glance on how TPNG patterns change before complications emerge, paving the way for future studies on prediction biomarkers and potentially disease mechanisms.

Performance of prediction models for nephropathy in people with type 2 diabetes: systematic review and external validation study.

OBJECTIVES: To identify and assess the quality and accuracy of prognostic models for nephropathy and to validate these models in external cohorts of people with type 2 diabetes. DESIGN: Systematic review and external validation. DATA SOURCES: PubMed and Embase. ELIGIBILITY CRITERIA: Studies describing the development of a model to predict the risk of nephropathy, applicable to people with type 2 diabetes. METHODS: Screening, data extraction, and risk of bias assessment were done in duplicate. Eligible models were externally validated in the Hoorn Diabetes Care System (DCS) cohort (n=11 450) for the same outcomes for which they were developed. Risks of nephropathy were calculated and compared with observed risk over 2, 5, and 10 years of follow-up. Model performance was assessed based on intercept adjusted calibration and discrimination (Harrell's C statistic). RESULTS: 41 studies included in the systematic review reported 64 models, 46 of which were developed in a population with diabetes and 18 in the general population including diabetes as a predictor. The predicted outcomes included albuminuria, diabetic kidney disease, chronic kidney disease (general population), and end stage renal disease. The reported apparent discrimination of the 46 models varied considerably across the different predicted outcomes, from 0.60 (95% confidence interval 0.56 to 0.64) to 0.99 (not available) for the models developed in a diabetes population and from 0.59 (not available) to 0.96 (0.95 to 0.97) for the models developed in the general population. Calibration was reported in 31 of the 41 studies, and the models were generally well calibrated. 21 of the 64 retrieved models were externally validated in the Hoorn DCS cohort for predicting risk of albuminuria, diabetic kidney disease, and chronic kidney disease, with considerable variation in performance across prediction horizons and models. For all three outcomes, however, at least two models had C statistics >0.8, indicating excellent discrimination. In a secondary external validation in GoDARTS (Genetics of Diabetes Audit and Research in Tayside Scotland), models developed for diabetic kidney disease outperformed those for chronic kidney disease. Models were generally well calibrated across all three prediction horizons. CONCLUSIONS: This study identified multiple prediction models to predict albuminuria, diabetic kidney disease, chronic kidney disease, and end stage renal disease. In the external validation, discrimination and calibration for albuminuria, diabetic kidney disease, and chronic kidney disease varied considerably across prediction horizons and models. For each outcome, however, specific models showed good discrimination and calibration across the three prediction horizons, with clinically accessible predictors, making them applicable in a clinical setting. SYSTEMATIC REVIEW REGISTRATION: PROSPERO CRD42020192831.

Long RNA Sequencing and Ribosome Profiling of Inflamed ß-Cells Reveal an Extensive Translatome Landscape.

Type 1 diabetes (T1D) is an autoimmune disease characterized by autoreactive T cell-mediated destruction of the insulin-producing pancreatic ß-cells. Increasing evidence suggest that the ß-cells themselves contribute to their own destruction by generating neoantigens through the production of aberrant or modified proteins that escape central tolerance. We recently demonstrated that ribosomal infidelity amplified by stress could lead to the generation of neoantigens in human ß-cells, emphasizing the participation of nonconventional translation events in autoimmunity, as occurring in cancer or virus-infected tissues. Using a transcriptome-wide profiling approach to map translation initiation start sites in human ß-cells under standard and inflammatory conditions, we identify a completely new set of polypeptides derived from noncanonical start sites and translation initiation within long noncoding RNA. Our data underline the extreme diversity of the ß-cell translatome and may reveal new functional biomarkers for ß-cell distress, disease prediction and progression, and therapeutic intervention in T1D.